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Within the framework of a (1+1) -dimensional model which mimics high energy QCD, we study the 
behavior of the cross sections for inclusive and diffractive deep inelastic 7*/i scattering cross sections. 
We analyze the cases of both fixed and running coupling within the mean field approximation, in 
which the evolution of the scattering amplitude is described by the Balitsky-Kovchegov equation, 
(jj ■ and also through the pomeron loop equations, which include in the evolution the gluon number 

' fluctuations. In the diffractive case, similarly to the inclusive one, the suppression of the diffusive 

scaling, as a consequence of the inclusion of the running of the coupling, is observed. 
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O ; I. INTRODUCTION 

It is well known that the high energy regime of the Quantum Chromodynamics (QCD) is described by non-linear 
43 ■ evolution equations [1-19]. At the level of scattering amplitudes, and in the framework of the dipole picture [20- 
22] , the most general ones are the so called pomeron loop equations [23-26] , which correspond to a generalization of 
the Balitsky-JIMWLK hierarchy [7-14, 16-19], by including the gluon number fluctuations. If one performs a mean 
field approximation, this infinite set of equations reduces to a single closed equation for the scattering amplitude 
of one dipole with a hadronic target, the Balitsky-Kovchegov (BK) equation [16, 27, 28], the simplest of the non- 
■ linear equations for the scattering amplitudes in QCD at high energy. This equation admits [29-31] travelling wave 
solutions, which have become a natural explanation for the geometric scaling — first observed in the HERA data for 
electron-proton deep inelastic scattering [32, 33] — and, being a mean field version of the complete hierarchy, neglects 
ON ' the effects of the fluctuations. At least in the fixed coupling case, from the correspondence between high energy 
QCD and reaction diffusion processes, one of the consequences of the gluon number fluctuations in the evolution of 
the dipole scattering amplitudes is, at very high energies, the replacement of the geometric scaling [32, 33], by the 
diffusive scaling [34]. 

Fluctuation effects have not been observed in the experimental data yet. Besides, the only few phenomenological 
studies have been inconclusive with respect to their presence in the current experiments [35-38]. Their physical 
consequences in the high-energy evolution in QCD for the phenomenology were first analyzed in Ref. [34], where their 
effects in the behaviour of inclusive and diffractive cross sections for deep inelastic lepton-hadron scattering (DIS) 
were studied. They found, for example, that, within the high energy regime, all the amplitudes or cross sections show 
diffusive scaling, that is, they depend upon the photon virtuality Q 2 and the total rapidity Y through the variable 
ln(Q 2 / (Q 2 )), where {Q 2 S ) is the (average) hadron saturation momentum. 

Our current knowledge on the consequences of the fluctuations comes only from the correspondence between high 
energy QCD and statistical physics; because of the complexity of the pomeron loop equations, the properties of the 
solutions are known only after some approximations, in asymptotic regimes and at fixed coupling [23]. On the other 
hand, in the last few years one observed an important progress in the inclusion of next-to-leading order (NLO) effects 
in the non-linear mean field BK equation. In particular, one can cite the explicit calculation of the running coupling 
effects [39-45] and its successful use in the description of HERA and RHIC data [46-48] . Unfortunately, because of the 
complexity of the pomeron loop equations, the inclusion of such NLO effects in these equations turn to be a very hard 
task. The difficulty of dealing with these equations, even in the fixed coupling case, inspired other ways of investigation 
of high energy evolution in QCD, in particular through particle models with a smaller number of dimensions [49-56]. 
Among them, the (l+l)-dimensional model presented in Ref. [54] has shown to mimic fixed impact parameter high 
energy QCD with fixed coupling constant. Its generalization to the case with the running coupling was done in Ref. 
[55]. In such version, the model could provide, for the first time, the study of both running coupling and fluctuations 
effects, taken into account simultaneously, in the high energy evolution of scattering amplitudes. The main conclusion 
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presented by the authors was the strong suppression of the pomeron loop (fluctuation) effects due to the running of 
the coupling, up to rapidity Y ~ 200, that is, well beyond the energies of interest for the phenomenology in QCD. 
The dynamics is similar to the respective prediction of the mean-field approximation with running coupling, the 
property of (approximate) geometric scaling being preserved for the average scattering amplitude. This result is in 
sharp contrast with the fixed coupling results, which show the emergence of the diffusive scaling. 

In this paper we present an investigation of the effects of both pomeron loops and running coupling, taken into 
account simultaneously, on the cross sections for inclusive 7* ft and, for the first time, on diffractive deep inelastic 
scattering (DIS), within the framework of the toy model [55]. In Sections II and III we present some important aspects 
of lcpton-hadron DIS, specifically an overview of kinematics and the description of the dipole picture of the inclusive 
and diffractive 7* ft, scattering. Section IV is devoted to an overview of the one-dimensional model in the running 
coupling case. In particular, we present the resulting evolution equations for the scattering amplitudes and the main 
features of their evolution. In Section V we present our results, with the study of the behaviour of the cross sections 
in both fixed and running coupling cases, and the conclusions are presented in Section VI. 



II. INCLUSIVE VIRTUAL PHOTON-HADRON DIS 



This process is described by the reaction l(k) + h(P) —¥ l(k')+X(P x ), where I refers to the lepton (with momentum 
k in the initial state and k' in the final one), ft to the incoming hadron (with momentum P) and X is the generic 
hadronic final state (with momentum Px)- Processes described by the reaction above are called inclusive, because 
only the lepton is measured in the final state. In the specific case where the lepton is an electron, its interaction with 
the hadron is mediated by a virtual photon with virtuality Q 2 = —q 2 = (k — k') 2 . If one looks at the -f*h — > X, in 
inclusive DIS all what is known from the final hadronic state X is that it has an invariant mass squared W 2 = (P + q) 2 , 
which is the center-of-mass energy of the 7* ft system. Another important definition is that of the Bjorken variable, or 
Bjorken-x, given by XBj = Q 2 / (Q 2 + IV 2 ); from it, one sees that, for fixed values of Q 2 , when one increases the energy 
IV 2 , XBj decreases and the high energy limit corresponds to the small-xej limit. The total rapidity of the process is 
defined as Y = ln(l/xBj)- 

At small-XBj, the 7* ft, process can be described in a convenient frame, the so-called dipole frame, in which the 
hadron carries most of the total energy, but the virtual photon has enough energy to split into a quark-antiquark (qq) 
pair, or a dipole. This dipole, then, interacts with the hadron. The dissociation of the virtual photon into the color 
dipole takes place long before the scattering, and the dipole evolves through soft gluon radiation until it meets the 
hadron (at the time of scattering) and scatters off the color fields therein. Exactly as it was done in [34], the present 
analysis will be restricted to the leading logarithm approximation, in which the evolution consists of the emission of 
soft gluons, carrying a small fraction x^j < 1 of the longitudinal momentum of their parent parton. In the limit 
N c — > 00, a gluon can be effectively replaced by a pointlike quark-antiquark pair in a color octet state, and a soft 
gluon emission from a color dipole can be described as the splitting of the original dipole into two new dipoles with a 
common leg. In this picture, the original qq pair produced by the dissociation of the virtual photon evolves through 
successive dipole splittings and becomes an oniura — i.e., a collection of dipoles — at the time of scattering. This is the 
Mueller's dipole picture [20-22]. 

Using the formalism developed in [34], one finds that the differential cross-section for onium-hadron scattering at 
fixed impact parameter is given by 

^(r,b,Y) = 2ReA(x,y;Y), (1) 

where A is the amplitude for the elastic scattering, b = (x + y)/2 and r are the impact parameter and the transverse 
size of the original dipole and x and y its transverse coordinates. 

In such high energy approximation, the DIS cross section for the inclusive virtual photon-hadron (7* ft) scattering 
can be expressed as 



<Ki5t (V,Q 2 )= / dv [d 2 r ]T \r a (r,v;Q 2 )\ 2 2Re A(x,y;Y), 

J ° J a=L.T 



(2) 



d 2 b 

where \ip^^ L \ 2 are the probability densities for the qq dissociation of a virtual photon with transversal (T) or longitu- 
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dinal (L) polarization, obtained from perturbative QED [3, 57], given by 

I*tM;Q 2 )| 2 = ^^E^i^ + C 1 "") 2 ] W(^) + "»SX?(<5,r)} (3) 

q 

\* L (r,v;Q 2 )\ 2 = ^^^{4(3V(l-^(0,r)}, (4) 

q 

where Q q = v(l — v)Q 2 + in 2 , m q is the mass of the quark with flavour q, Ko,i are the Mc Donald functions of rank 
zero and one, respectively, and v is the fraction of the photon longitudinal momentum carried by the quark. 

Expression (2) is a priori frame- independent, but, the inclusive cross section is most simply evaluated in the frame 
where almost all the total rapidity Y is carried by the hadron (the target) and the projectile is an elementary dipole. 
In this case, A(x,y;Y) = (T(x,y)) Y [34] and 



(Y,Q 2 ) = f&vftfr J2 \iPUv,r;Q)\ 2 2Rc(T(x,y)) Y , (5) 



where (T(x,y)) Y is the (average) one dipole-hadron scattering amplitude, the brackets meaning the average over 
the target configurations. Here, we are interested in the high-energy limit of the DIS cross sections at fixed im- 
pact parameter. We assume that the dependence on b can be factorized into a profile function S(b), according to 
(T(x, y)) Y = S(b) (T(r)) Y (r = \x — y\ is the dipole size), where the integral (To = J dbS(b) would provide an overall 
normalization factor of order of the transverse area of proton. Since the dependence on b in such an approximation 
results completely decoupled, in the following we simply set S(b) — 1, assuming the integration over b extended up 
to 6 max providing the correct normalization of the cross section. The differential inclusive cross section reads 



pl poo 

(Y,Q 2 )=4n dv dr V \^ a (v,r; Q)\ 2 r Re (T(r)) Y 
d b Jo Jo a=L T 



(6) 



As it will be convenient for our purposes, we can write (T(r)) Y = (T(x)) Y = (T X ) Y , where x = ln(l/r 2 Q§) represents 
r in logarithmic units 1 (Qo is a scale of reference introduced by the initial conditions at low energy). The total cross 
section, then, takes the form 



d(J tot 

d 2 b 



(^Q 2 ) = 7^/ d W A *r~ x E \r a (v,x;Q)\ 2 Re(T x ) Y . 

WO JO J-OO n _J rp 



(7) 

a— L,T 



III. DIFFRACTIVE DIS 

Part of the DIS events are diffractive. In such events, described by the reaction 7*/i4 Xh, the final states contain 
an intact scattered hadron h and a diffractive hadronic state X separated by a rapidity gap Y gap = ln(l/a;p), where 
xp = XBj/P and j3 is related to the diffractive invariant mass Mx by j3 = Q 2 / (Q 2 + M x ). It is straightforward to see 
that the difference between the total rapidity Y and the rapidity gap Y gap is Y — Y gap = ln(l//?). 

The cross section for the diffractive process reads (Y gap denotes the minimal rapidity gap) 



cr 7 f 1 f 

^(Y,Y gap ,Q 2 ) = dvjd 2 r Y mv,r;Q)\ 2 P dis (b,r;Y,Y gap ). (8) 



Since the whole process can be factorized, for our purposes it will be enough to ignore the electomagnetic process 
(the splitting the virtual photon into the qq dipole) and focus only on the onium-hadron (Oh) scattering. More 
specifically, we will be interested in the quantity Pdiff , which is the probability for diffractive onium-hadron scattering 
(Oh — > Xh), and corresponds to the differential cross-section for onium-hadron scattering at fixed impact parameter: 

^ (r, b, Y, Y gap ) = P diS (x, y; Y, Y gap ). (9) 



1 The variable x should not be confused with the bold-faced x, which represents a vector in the transverse plane in the picture of DIS. 



4 



An explicit formula for this probability has been obtained within the lightconc wavefunction formalism in [34] , in a 
special frame, in which Y gap coincides with the rapidity Y of the target hadron. This choice of the frame is important 
because it avoids one to deal explicitly with final state interactions. The resulting formula is given by 

P d i ff (x, y; Y, Y gap ) -+ P diff (x, y; Y, Y ) = ]T P({N}; Y - Y ) 

{N} 

which has the following meaning: starting from an original dipole (x, y), after an evolution Y — Y there is a probability 
density P({N};Y — Y ) for a given configuration of N dipoles to be produced. S Xi = 1 — T Xi is the S-matrix for 
the scattering between the ith dipole (with logarithmic size Xi) in the projectile and a given configuration of the 
target, T Xi being the corresponding T-matrix. The symbol X){at} represents the sum over all the configurations of 
the projectile with N dipoles. Again, the notation {-) Y denotes the average over the ensemble of color fields in the 
target. 

Our main aim is to investigate, for the first time, the behavior of the diffractive probability (10) with increasing 
rapidity Y in the presence of fluctuations and running coupling effects. This requires the description of the rapidity 
evolution of the dipole-hadron scattering amplitude T x , as well as the probability density P({N}; Y — Yq) by taking 
into account both effects simultaneously, which is still a prohibitive task in full QCD. However, a convenient way of 
doing this is through the model presented in [55], whose main features we will briefly describe below. 



IP* 



(10) 



IV. (l+l)-DIMENSIONAL MODEL FOR HIGH ENERGY QCD 

The toy model [55] is a (l+l)-dimensional stochastic particle model, where one of the dimensions refers to the 
total rapidity separation Y between two hadronic systems which undergo evolution and scattering (and plays the 
role of time in the evolution), while the other one (the spatial dimension) is the position of the particle along an 
infinite one-dimensional axis, the x-axis, which, in analogy with the dipole picture of QCD [20, 22], corresponds to 
the logarithm of the inverse size of a dipole, as defined in Section II. 



A. The structure of model 



In this toy model, a system of particles (which corresponds to a given hadronic system) is specified by their 
distributions along the one-dimensional x-axis. In order to describe a scattering problem, one considers two such 
systems (projectile and target) which scatter off each other along a given collision axis (which is transverse to the 
x-axis) and assumes that each particle of the projectile can scatter elastically with any particle of the target. The total 
rapidity Y of the process is divided between the right mover system (R), the projectile, which has rapidity SY = Y — Y 
and the left mover system (L), the target, which has rapidity —Y . Let P R [n(x R ),Y — Y ] and Pi\m{xi), Yq\ be the 
probability densities to find given configurations in the two systems, these being described as functions of the densities 
of particles at the point x. The average S-matrix is given by 

(S) Y = J VnVmP R [n(x R ),Y -Yo}P L [m(x L ) 7 Y }S[n(x R ),m(x L )}. (11) 

Here, S[n(x R ) , m(x^)] is the S'-matrix associated with a given pair of configurations and the (• • •) symbol represents 
the average over all possible configurations {n(x R )}, {to(:el)}- This 'event- by- event' S'-matrix is given by 



S[n, to] = exp 



J dx R dx L n(x R )m(x L ) lna(x R \x L ) 



(12) 



where <j{x r \xl) = 1 — t(xr|xl) is the S-matrix for the scattering of two elementary particles of logarithmic sizes xr 
and xl, and t(xr|xl) the corresponding T-matrix (0 < t(xr|xl) < 1). 

The probability densities obey the following evolution equation (the details of the evolution can be found in refer- 
ences [54, 55]) 



dP[n(x),Y] 
d~Y 



J dzf z [n(x)-S(x-z)]P[n(x)-S(x-z),Y}- J dz f z [n(x)] P[n(x), Y], (13) 



2 Here we follow the same notation used in [55] 
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where f z [n(x)] is the probability per unit rapidity to find an extra particle with logarithmic size z after an evolution 
step (after a small increment in rapidity, only one extra particle can be emitted), given that the initial configuration 
of the evolved system was n{x). The functional form of the "deposit" rate density f z [n(x)] can be found by assuming 
Lorentz invariance, and one gets 



fzHx)] = 



T z [n(x) 
a(z) 



(14) 



where T z [n(x)] is the T-matrix for the scattering of a particle of logarithm size z off a system with a given configuration 
n(x), and is given by 



T z [n(x)] = 1-exp 



/ 



dxn(x) Inertia:) 



(15) 



and a(z) is the coupling parameter. In the case of running coupling, a(z) = l/f3z — in such a way to mimic the one- 
loop running coupling of QCD — , with j3 being the analog of the one-loop beta function of QCD. Another important 
feature of the model is the specification of the explicit form of the elementary particle-particle scattering amplitude 
t(x\u), which, in analogy with the corresponding quantity in QCD (the amplitude for dipolc-dipole scattering), is 
chosen as 



r(x\y) = a(x)a(y) exp(-|a; - y\) = a(x)a(y)K(x, y) = a x a y K xy . 



(16) 



With the above expressions at hand, one can now present the evolution equations for any observable. In particular, 
since we want to describe the cross section (7), we will present the resulting equations for the scattering amplitudes. 



B. Evolution of the amplitudes 

Let us consider a generic observable O which depends on the configuration of the particles in the system. If one 
evaluates its average value at rapidity Y, one gets a measurable quantity, given by 

(0) Y = J VnP[n(x),Y}0[n(x)}, (17) 

where here we mean that the average is taken over all the configurations of the right mover, that is, P[n(x),Y] = 
Pji[n{x),Y] (the left mover consists in a given configuration of particles). 

By using Eqs.(13) and (17), it is straightforward to obtain the evolution equation for any physical observable: 



9 (Oh 
dY 



= Jdz (f z [n(x)} {0[n(x) - 5(x - z)] - 0[n(x)}}) Y . (18) 



If the observable is the amplitude for the scattering between a projectile which consists of a single particle of a 
given logarithmic size x and a generic target, one has (the average over Y is implicit) 



d(T x ) 
dY 



a x j ' dzK xz (T z (l-T x )), (19) 



which is not a closed equation for (T x ), but the first equation of an infinite hierarchy. This equation is analogous 
to the first equation of the Balitsky-JIMWLK hierarchy, which is identical to the first of the pomeron loop (PL) 
equations (extended to running coupling): it has a linear term, proportional to (T), and a non-linear (quadratic) 
term, proportional to T 2 . The term corresponding to the particle number fluctuations appears only in the second 
equation of the hierarchy, which reads 



d(T x T y ) 



a x f dzK xz (T z T y (l - T x )) + a y J dzK yz (T Z T X (1 - T y )) 



dY 

+ a x a y J ' dza z K xz K yz (T z {l-T x ){l-T y )) , (20) 
the fluctuation term being the one proportional to (T) in the second line of the above equation. 
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In the mean field approximation (MFA) , the whole hierarchy reduces to a single closed equation, which is obtained 
by making (TT) = (T) (T) in Eq.(19), resulting in the analogous to the (running coupling) Balitsky-Kovchegov (BK) 
equation 

^ =«x/ K xz [(T z ) - (T z ) (T x )} . (21) 

z 

The evolution equations with fixed coupling can be straightforwardly obtained by simply making a = constant. Now, 
we will make a brief review of the main aspects of the evolution of the average amplitude (T), in both fixed and 
running coupling cases. 

(i) Fixed coupling case [54]: 

1. From the similarity with the BK equation, Eq.(21) admits the so called travelling wave solutions, which means 
that, at very large values of rapidity, the amplitude depends on x and Y through the scaling variable x — x s (Y), 
that is, the amplitude T is a front which interpolates between 1 and and, as Y increases this front gets simply 
translated towards larger values of x, without being distorted. The function x s (Y) is the saturation scale, which 
naturally emerges from the non-linear evolution; it separates between the dense target region, x x s , where 
T = 1, and the dilute target region, x ^ x s , where T decreases exponentially. It is also an increasing function 
of rapidity, the analogous to the (logarithm of the) saturation momentum which emerges from the non-linear 
evolution in QCD, \n(Q 2 /Q^). It can be also defined as the position of the front, that is, the line along which 
the amplitude is constant and of 0(1) (it is usual to choose T(x = x s (Y),Y) = 1/2). The dependence on the 
combined variable x — x s is the so-called geometric scaling [32, 33], which is valid in a window which grows with 
increasing rapidity like oc Y" 1 / 2 . 

2. With the inclusion of the fluctuation effects, the differences with respect to the mean field analysis are very 
significant. From a given initial condition (at Y = 0), the evolution up to Y generates a statistical ensemble 
of fronts, which have the same form, but differ from each other by their respective front positions x s , and this 
position is now a random variable. In the fixed coupling case, to a very good approximation, the distribution 
of x s with Y is a Gaussian, with both the expectation value (x s ) and the dispersion a 2 rising linearly with Y. 
The individual fronts exhibit geometric scaling, but only over a compact region, in contrast with the mean field 
amplitude, for which the scaling window is ever increasing with Y. The average amplitude (T) is obtained by 
averaging over the ensemble, and one gets that the geometric scaling property of the individual fronts is washed 
out by the dispersion of the fronts, and is replaced, at sufficiently large Y, by the diffusive scaling. 

(ii) Running coupling case [55]: 

1. In the MFA, geometric scaling is also present at asymptotic rapidities in the evolution of the amplitude, but the 
front formation is delayed in comparison with the fixed coupling case: the window for geometric scaling grows 
with increasing rapidity like oc Y 1 / 6 . 

2. After the inclusion of the fluctuations, the growth of the dispersion with Y is suppressed, and one has a 2 oc y/Y. 
Besides, the influence of the fluctuations is strongly suppressed, remaining negligible for all the rapidities of 
interest. In particular, the average amplitude exhibits approximate geometric scaling: the deviation from 
geometric scaling with increasing Y is too small. 



V. RESULTS 



Now we can study the consequences of the properties of the scattering amplitudes discussed above on the behavior 
of the cross section (7) and the diffractive probability (10) with increasing Y. Concerning the inclusive cross section 
(7), our aim here is only to reproduce the results obtained in [34] and [55], respectively, at fixed and running coupling. 
The diffractive case is our main result: in the high energy limit, it has been shown that, in the fixed coupling case, 
diffractive cross section exhibits diffusive scaling [34] . Here, for the first time, we study the behavior of this quantity 
with increasing energy in the presence of both fluctuation and running coupling effects. 

For the purposes cited above, we must use as the input for the average scattering amplitude for particle-hadron 
(dipole- hadron) the solution of Eqs.(19) (when fluctuations are included) and (21) (when fluctuations are not included, 
that is, in the MFA). The parameters which enter into the expression for the cross section (7) must be fixed: we set 
Qq = 1 GeV 2 , the electromagnetic coupling constant a e = 1/137, the number of colors N c = 3, and only light quarks 
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(u, d, s) enter into this analysis, with zero masses. In both inclusive and diffractive cases, we perform the analysis 
first considering the fixed coupling case, in both mean field appoximation (MFA) and with fluctuations included and 
then we generalize it by doing the same in the case with running coupling. In the specific case of diffractive onium- 
hadron scattering, in the evaluation of the probabilities P({N}; 6Y) for a given configuration of the projectile onium 
at rapidity 5Y, as well as the averages over all the target configurations, we follow the same procedure described in 
[54]. 

A. Fixed coupling case 

The results in the case of fixed coupling (FC) are shown in Figures 1 and 2. The value of the coupling constant 
is chosen to be a = 0.2. Fig. 1 presents the DIS inclusive cross section as a function of the variable Q 2 / (Qf), for 
different values of rapidity, up to Y — 100. One should remember that the average saturation momentum, (Q^(Y)), 
is related to the average saturation scale (x s (Y)): (x s (Y)) = \n({Q 2 s (Y)) = ln(l/ (r s ) 2 Ql). In the MFA (left 
plot), one clearly sees the geometric scaling, as well as the growth of its window as rapidity increases from Y = 0. For 
values of rapidity values Y ?Z 30, the curves for the cross section have the same shape and they depend only on the 
scaling variable Q 2 / (Q 2 )- After the inclusion of the fluctuations (right plot), the curves deviate from the mean field 
behaviour (and thus from geometric scaling) as Y increases. These FC results for the inclusive cross section reflect the 
corresponding behaviour of the scattering amplitude and are consistent with the ones already obtained in the QCD 
framework [34]. 




Q 2 /(Qs) 2 Q 2 /(Q S ) 2 

FIG. 1: Fixed coupling results for various rapidities as a function of Q 2 /(Q 2 ), in the MFA (left plot) and with the inclusion of 
the fluctuations (right plot), up to rapidity Y = 100. 

In Fig, 2 the diffractive probability for onium-hadron DIS is shown as a function of the geometric scaling variable 
(r s ) 2 /r 2 = e~^ x '^~ x ^ for different values of the total rapidity interval Y. The rapidity interval of the projectile onium, 
SY = Y — Yo, is kept fixed at a small value (5Y — 1), to ensure that the projectile is a dilute system, consisting 
of a small number of particles (dipoles). In the MFA, geometric scaling is reached at very large values of Y. When 
fluctuations are included, one observes that, similarly to the inclusive cross section, geometric scaling breaks down 
and Pdif f exhibits diffusive scaling. This result is consistent with what has been found in [34] 

B. Running coupling case 

Now, we proceed with a generalization of the previous case, by taking into account the running of the coupling, 
given by a x = l/f3x, with (3 chosen to be 0.72. The results are shown in Figure 3, where the cross section (7) is 
represented as a function of the variable Q 2 / (Ql), for different values of rapidity, up to Y — 200. In the MFA, one 
can observe geometric scaling, but it is reached at larger values of rapidity in comparison with the FC case. This 
reflects the corresponding behaviour of the scattering amplitude, for which the formation of the front in the RC case 
is delayed. With the inclusion of the fluctuations, one can observe that the increasing dispersion present in the FC 
case is strongly suppressed and one has an approximate geometric scaling, since the different curves have quite small 
deviations from each other when increasing rapidity, resulting in a behavior very similar to the MFA (with running 
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(r s ) 2 /r 2 (r s ) 2 /r 2 

FIG. 2: Diffractive probability for onium-hadron scattering as a function of the scaling variable (r 2 ) /r 2 : fixed coupling results 
for various rapidities, in the MFA (left plot) and with fluctuations (right plot), up to rapidity Y = 95. 

coupling) . Then, the high energy behavior of the inclusive cross section reflects the corresponding behaviour of the 
average particle (dipole) scattering amplitude in the running coupling expected. 




1 10 2 10 4 1 10 2 10 4 

Q 2 /(Qs} 2 Q 2 /(Q S } 2 

FIG. 3: Running coupling results for various rapidities as a function of Q 2 /(Q 2 }, in the MFA (left plot) and with the inclusion 
of the fluctuations (right plot), up to rapidity Y = 200. 

Our next step is to investigate if the suppression obtained in the inclusive case, due to both fluctuation and running 
coupling effects, holds also for the diffractive probability Pdiff(r, Y, SY). This answer is the main result of this paper. 
First, from the left plot in Fig. 4 we can see that, in the MFA, geometric scaling is observed, as expected, but is 
reached faster than in the FC case, at smaller values of rapidity (now SY = 2). Finally, in the right plot we present, 
for the first time, the study of the behavior of the diffractive probability in the presence of both fluctuation and RC 
effects. The suppression of fluctuations exists and is as strong as in the MFA case. Therefore, in diffractive DIS, 
within the framework of the toy model for high energy QCD, fluctuations are strongly suppressed by the running of 
the coupling and diffusive scaling of the cross sections, predicted in the FC case, is washed out. 

VI. CONCLUSIONS 

In this paper we have investigated the high energy behavior of the total cross section for virtual-photon-hadron 
DIS and for onium-hadron diffractive DIS within the framework of the (l+l)-dimensional model [55], which provides 
a way to study, at fixed impact parameter, the effects of the particle number fluctuations and running coupling, 
taken into account simultaneously. In the fixed coupling case, the results are consistent with those obtained in the 
framework of QCD [34], that is, the geometric scaling which is present in the mean field approximation at large values 
of rapidity, is completely washed out when fluctuations are taken into account. 
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FIG. 4: Diffractive probability for onium-hadron scattering as a function of the scaling variable (r 2 ) /r 2 : running coupling 
results for various rapidities, in the MFA (left plot) and with fluctuations (right plot), up to rapidity Y = 190. 



By generalizing the analysis done in [34], through the inclusion of running coupling effects, we have reproduced the 
results obtained in [55] for the inclusive virtual photon-hadron cross section: the behaviors of this cross section with 
and without fluctuations are similar, this observable presenting approximate geometric scaling, which means that the 
running of the coupling suppresses the fluctuation effects at asymptotic rapidities. In the diffractive onium-hadron 
scattering, the diffractive probability exhibits geometric scaling in the MFA. When fluctuations are included, diffusive 
scaling is seen in the fixed coupling case, while that, in the running coupling case, geometric scaling is present and 
reached at smaller values of rapidity Y than in the case without fluctuations. 

This suggests that the mean field treatment with running coupling would be enough to study not only the inclusive 
lepton-hadron DIS, but also the diffractive DIS, for all the energies available at present and to be available in a near 
future. The toy model also allows the investigation of the other processes which, in the framework of QCD, admit a 
dipole factorization. Thus, it would be interesting to apply it to such processes, in particular less inclusive ones, in 
order to investigate if the suppression of fluctuations by running coupling effects remains present. 
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